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ABSTRACT 

lemorandum  presents  a  usage  manual  for  an  IBM  7090  computer  program 
which  uses  a  Monte  Carlo  simulation  to  determine  the  probability  of  destroying  a 
rectangular  target  with  air-to  ground  gunnery.  The  effect  of  correlation  of  suc¬ 
cessive  aim  points  is  considered.  It  is  assumed  that  the  aim  point  is  distributed 
with  a  bivariate  normal  aiming  error  and  that  the  individual  rounds  are  distributed 
with  an  independent  bivariate  normal  ballistic  dispersion.  A  flow  chart,  a  listing 
of  the  FORTRAN  program  and  a  sample  problem  are  included. 
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I.  INTRODUCTION 

Hitherto,  models  for  air-to-ground  gunnery  have  been  of  two  kinds.  The  first 
assumes  that  the  aim  points  for  successive  rounds  are  uncorrelated  and  independent, 
i.e.  that  the  correlation  of  successive  aim  points  is  zero.  The  second  model  is  the 
salvo  type  which  assumes  that  successive  aim  points  are  completely  correlated.  It 
is  well  known  that  both  of  these  models  are  unsatisfactory  desc^ij  Tons  of  air-to- 
ground  gunnery  since  the  aim  point  does  wander  and  the  correlation  is  neither  Onor 

I,  but  some  intermediate  value.  These  two  models  have  been  used  because  the 
more  general  case  defies  simple  mathematical  analyses.  However,  with  the  intro¬ 
duction  of  high  speed  computers  and  Monte  Carlo  techniques,  it  is  possible  to 
simulate  the  more  general  case. 

II.  GENERAL  DESCRIPTION 

This  program  is  a  Monte  Carlo  simulation  for  use  on  the  IBM  7090  computer. 

It  determines  the  probability  of  destroying  a  rectangular  target  with  air-to-ground 
gunnery  on  the  basis  of  the  following  assumptions: 

a.  The  aiming  error  and  the  ballistic  dispersion  are  bivariate  normal  and 
are  independent  in  the  range  coordinate  (along  the  flight  of  path  of  the  aircraft) 
and  the  deflection  coordinate  (normal  to  the  flight  path  in  the  horizontal  plane). 

b.  There  is  a  conditional  kill  probability  (probability  of  kill  given  a  hit)  as¬ 
sociated  with  each  round,  and  there  is  no  cumulative  damage  effect. 

c.  The  correlation  coefficient  between  the  n^and  m^1  round  is  a  m  ^  in 

In-m  I 

the  range  coordinate  and  bn  in  the  deflection  coordinate  where  a  and  b 
represent  the  correlation  constants  (Note:  when  a  and  b  are  zero  or  one,  the 
process  degenerates  into  the  independent  or  salvo  cases,  respectively). 

d.  The  mean  aim  point  is  at  the  center  of  the  target  which  is  rectangular 
and  is  located  with  one  axis  parallel  to  the  flight  path. 

In  addition  to  the  parameters  specifying  the  target  size,  aiming  and  bal¬ 
listic  dispersion,  correlation  and  conditional  kill  probability,  the  program  inputs 
include  the  aircraft  speed,  the  firing  rate  of  the  gun,  the  slant  range  from  the 
target  at  the  commencement  of  the  firing  run,  the  maximum  number  of  rounds 
per  pass,  the  number  of  Monte  Carlo  runs  to  be  made  and  the  Increment  in  the 
number  of  rounds  at  which  target  kill  prolabilities  should  be  determined. 

Ill  MATHEMATICAL  DERIVATION 

Let  x^  and  y  be  random  variables  denoting  the  position  of  the  aim  ooint 
th 

of  the  n  ‘  round  in  the  plane  normal  to  the  line  of  sight  between  the  aircraft  and 
the  target.  Let  a  and  be  independent  (Russian  random  variables  with 

standard  deviation  1  and  mean  0.  Let 
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IV.  METHOD  OF  SOLUTION 


The  inputs  for  the  program  consist  of  the  following  parameters: 

Address* 

Symbol 

Description 

1 

a 

Correlation  constant  in  the  range  coordinate  (Osasl) 

2 

b 

Correlation  constant  in  the  deflection  coordinate, 
(0sb*l) 

3 

°R 

Standard  deviation  of  aim  error  in  the  range  coordinate 
(mils) 

4 

°D 

Standard  deviation  of  aiiu  error  in  the  deflection  co¬ 
ordinate  (mils) 

5 

0R 

Standard  deviation  of  ballistic  dispersion  in  the  range 
coordinate  (mils) 

6 

0D 

Standard  deviation  of  ballistic  dispersion  in  the  de¬ 
flection  coordinate  (mils) 

7 

S 

Slant  range  at  commencement  of  firing  run  (ft.) 

8 

R 

Firing  rate  (rounds  per  minutey 

9 

c 

Aircraft  speed  (knots) 

10 

N 

Maximum  number  of  rounds  per  pass 

11 

P 

Conditional  kill  probability  of  a  round 

12 

L 

Target  length  (normal  to  flight  path  in  range  coordinate) 
(ft.) 

13 

W 

Target  width  (normal  to  flight  path  in  deflection  co¬ 
ordinate  )  (ft.) 

14 

F 

Number  of  Monte  Carlo  iterations 

15 

N 

r 

Number  of  dummy  passes  through  random  number 
generators 

16 

AN 

Increment  in  the  number  of  rounds  at  which  target  kill 
probabilities  are  determined. 

•Tile  parameter  addresses  are  explained  in  appendix  c,  Data  Subroutine. 


Unless  otherwise  specified,  the  following  parameters  are  assigned  the 


designated  values. 

Address  Parameter  Value 


1  a  0 

2  b  0 

9  c  0 

11  P  10 

14  F  1000 

16  AN  1 

15  Nr  0 
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The  program  flow  chart  and  FORTRAN  statements  are  shown  in  appendixes  A 
and  B.  ^he  subroutines  are  described  in  appendixes  C,  D  and  E.  The  solution  is 
obtained  in  the  following  manner: 

1.  Before  starting  the  first  iteration  for  the  first  data  set,  nr  dummy  passes 
are  made  through  the  random  number  generator. 


2.  For  each  Monte  Carlo  iteration: 


a.  Given  the  standard  deviations  of  aim  error  (o^  ,  c^)  and  the  aim 

st 

point  of  the  (n-1)  round,  two  Gaussian-distributed  (mean  =  0, 
standard  deviation  =  1)  random  numbers  (a,  $  )  are  selected,  and 

and  the  aim  point  of  the  nu  round  (R  ,  D  )  is  determined: 

n  n 

R  =  aR  .  +cD  (l-a2)i/2o  (1) 

n  n-1  R ' 

D  =  bD  ,  +ari(l-b2)1/2  8  (2) 

n  n-1  D 

Fc.  the  first  round  (n  =  1): 


R1  =  Or  a  (3) 

<4) 


b.  Given  the  target  dimensions  (L,  W),  the  firing  rate  (R),  the  air¬ 
craft  speed  (c)  and  the  initial  slant  range  (S),  the  half -target 
size  (in  mils)  is  determined  for  the  time  at  which  the  ntn  round 
is  fired: 


1/2  L 

n 


500L 


S  -  (n-l)c(~)(1.688) 


(5) 


1/2W 

n 


500W  _ 

S  -  (n-l)c(^)(1.688) 


(6) 


e.  The  previously  computed  aim  point  is  checked  to  det^-mine  if  it 
is  within  three  standard  deviations  of  ballistic  dispersion 
(3  r  ,  6  g)  from  the  target: 


|Rnl«.l/2Ln+3eR? 

(7) 

|Dj  s  1/2W^+30d  ? 

(8) 
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If  the  aim  point  is  not  within  three  standard  deviations  of  ballistic 
dispersion  from  the  target  in  either  coordinate,  it  is  assumed  that 

St 

the  round  misses  the  target  and  the  aim  point  for  the  (n  +  1)  round 
is  computed  (part  a, ). 

If  the  aim  point  is  within  three  standard  deviations  from  the  target 
in  both  coordinates,  a  Gaussian-distributed  random  number  (6)  is 
selected,  and  the  impact  point  of  the  round  in  the  range  coordinate  is 
determined  and  checked  to  ascertain  whether  the  round  falls  within 
the  limits  of  the  target  in  the  range  coordinates. 

lRn+/3R6i  <1/2Ln?  (9) 


If  the  round  does  not  fall  within  the  target  limits,  the  aim  point  for 
St 

the  (n  +  1)  round  is  computed  (parr  a.). 

If  the  nth  round  does  fall  within  the  target  limits  in  the  range  co¬ 
ordinate,  another  Gaussian-distributed  random  number  (e)  is 
selected  and  the  impact  point  of  the  round  in  the  deflection  co¬ 
ordinate  is  determined.  Another  check  is  then  made  to  determine 
if  the  round  hits  die  target: 

lDn  +  0De|<1/2W„?  <10> 

If  the  round  does  not  hit  the  target,  the  aim  point  for  the  fn  +  l)st 
round  is  determined  (part  a.). 

th 

If  the  n  round  does  hit  the  target,  a  uniformly  distributed  (between 
0  and  1)  random  number  (PP)  is  selected  and  compared  with  the 
conditional  kill  probability  to  determine  if  the  hit  results  in  a  target 
kill,  if  P  <  PP,  the  weapon  does  not  destroy  the  target  and  the  aim 
St 

point  of  the  (n  +  1)  weapon  is  computed  (part  a.}, 
th 

If  P  -  PP,  the  n  weapon  does  kill  the  target.  The  count  (N^)  of  the 

Monte  Carlo  iterations  for  which  the  target  is  destroyed  is  increased 
by  one.  In  addition,  a  counter  (JJ^)  for  the  next  highest  multiple  (i) 

of  the  increment  (dN)  in  the  number  of  rounds  for  which  the  probability 
is  to  be  determined  is  increased  by  one: 


Nh  +  1~Nh 

JJi  +  1  -Hi 


(11) 


When  a  round  has  resulted  in  a  kill  or  when  the  maximum  number  of 
rounds  (N)  has  been  fired  without  killing  the  target  the  entire  process 
is  repeated  until  F  Monte  Carlo  iterations  have  beer  completed. 
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3  At  conclusion,  the  estimated  probability  of  destroying  the  target  with  N 
rounds  is  determined: 

P(N)  =  Nh/F.  (12) 

The  probability  of  destroying  the  target  with  j  rounds  is  determined: 

n 

p(j)  =  .I1  JVF  (13) 

j=niN  (n  =  1,2,  .  ..N/AN)  (14) 


V.  USER’S  INSTRUCTIONS 

Input  flexibility  has  been  attained  by  allowing  the  user  to  vary  any  or  all  of 
the  parameters  in  a  computer  run.  There  is  no  programmed  limit  to  the  number 
of  data  sets  which  a  user  may  submit  in  a  run.  The  only  restriction  is  that  each 
data  set  must  terminate  with  one  blank  card,  and  the  last  set  in  the  run  must 
terminate  with  two  blank  cards. 

For  each  data  set  after  the  first,  the  user  need  submit  only  those  param¬ 
eter  values  in  a  set  that  are  different  from  those  in  the  previous  set.  This  is 
accomplished  by  identifying  each  input  parameter  by  its  address. 

VI  SAMPLE  PROBLEM 

A  target  whose  dimensions  in  the  plane  normal  to  the  flight  path  are 
10  feet  x  10  feet  is  to  be  attacked  by  a  500  KTAS  aircraft  with  a  firing  rate  of 
1000  rounds  per  minute.  The  aircraft  commences  firing  at  a  slant  range  of 
5000  feet  and  fires  a  single  burst  of  100  rounds  at  the  target.  The  probability 
of  kill  given  a  hit  is  0.2.  The  standard  deviations  of  aiming  error  are  6  and  4 
mils  in  range  and  deflection,  respectively,  and  the  standard  deviations  of  bal¬ 
listic  dispersion  are  2  and  1  mils  in  range  and  deflection.  The  correlation 
constants  in  both  range  and  deflection  are  assumed  to  be  0.  8.  What  is  the 
probability  of  killing  the  target  with  the  entire  burst  and  with  smaller  bursts 
in  increments  of  10  rounds?  A  Monte  Carlo  simulation  of  2000  iterations  is 
selected  The  input  statements  are  as  follows: 

Address  Value 

1  0.8 

2  0.8 

3  6 

4  4 

5  2 

6  1 

7  5000 

8  1000 
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(Continued  from  preceding  page.) 


Address 

9 

10 

11 

12 

13 

14 

15 

16 


Value 

500 

100 

0.2 

10 

10 

2000 

0 

10 


The  submittal  form  and  output  for  this  sample  problem  are  shown  in 
appendix  F . 

VII.  KEYPUNCH  INSTRUCTIONS 

Two  methods  of  punching  the  input  parameters  can  be  used.  In  the  first 
method,  one  card  is  used  for  each  parameter.  The  address  of  the  parameter 
is  punched  first,  is  followed  by  a  space  and  then  the  parameter  value  is 
punched  In  the  second  method,  the  packed  card  method,  several  parameters 
and  their  addresses  can  be  punched  on  a  single  card  (as  many  as  space  permits). 
In  each  case,  the  address  is  followed  by  its  parameter  value  which  is  in  turn 
followed  by  the  next  address  and  value.  Each  address  and  parameter  value 
must  be  separated  by  a  space  and  must  also  be  separated  from  the  value  pre¬ 
ceding  or  the  address  following  it. 

VIII.  OPERATOR'S  INSTRUCTIONS 

R.un  under  control  of  the  Bell  System  on  the  IBM  7090.  No  special  in¬ 
structions  are  required. 

IX.  TIMING 

This  program  requires  approximately  6  NF  milliseconds  running  time 
per  data  set.  where  N  is  the  number  of  rounds  per  pass  and  F  is  the  number 
of  Monte  Carlo  iterations. 
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APPENDIX  A 
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APPENDIX  B 

FORTRAN  STATEMENTS 


DIMENSION  0(165 >  JJ( 2000) ,FF < 2000 ) 

_ EOLU  l  VALENCE.  LDlJLUAl^l  0(2).B).(D(3).SIGR).(D(45.SIGD).ID(  5) .BETAR 

1 ) , (D( 6) ,BETAD) , (D( 7) ,S ) , ( D( B ) ,R) * ( D(9 ) ,C) , ( D( 10) , FN) « ( D( 11 )*P ) » 

A-0.0  ~  ~ 


_ _  B-Q.O _ 

F-1000.0 

_ Pj?-L».Q _ 

DN-1.0 

C-0.0  _ 

10  CALL  DATA  (O.INO) 

_ IF  1 1ND)  40,40 j20 

20  PRINT  30 
30  FORMAT  UH11 
CALL  ENDJOB 
40  PRINT  30 
II-DI 15) 


DO  45  1-1,11 _ _  _ __ _ 

CALL  RANUMB  (DUMMY) 

45  CALL  GRNUMB  (DUMMY) _ 

FLR2-500.0-FL/S 

WR2-500.C-W/S  _ 

V-i.688*C 

_  VR-60.0«V/(S*R)  _  _ _ 

N-FN/ON 

DO  50  1-1, N _ 

50  JJID-0 

__  NH-O _ _ 

JJJ-FN 

_ IJL-F _  _ _ _ 

DO  170  1*1,11 

CALL  GRNUMB  (ALPHA) _ _ 

CALL  GRNUMB  (GAMMA) 

RC-SIGR*ALPHA 

DC-SIGD-GAMMA 

_ DO  ISO  J-1,JJJ  _ _ _ 

60  IF  (VR)  80,70,80 

70  WN2-WR2 _ 

FLN2-FLR2 

GO  TO  90 _  _ _ _ 

80  WN2-WR2/ ( i.O-FLQATF ( J-l ) *VR > 
f LN2-FLR2/( l.O-FLOATF ( J-l )*VR ) 

90  i F ( AB  SF ( RC )— FLN2— 3 . 0-BETAR)  100,140,140 
100  I F  t  ABSF  (  DO-WN2-3 «0»8F TADS  110,140,140 
110  CALL  GRNUMB  (DELTA) 

IF _ (  ABSF((U>DELTA»BETAR)-FLN2 )  120, 140,140 

120  CALL  GRNUMB  (EPS) 

_  IF  ( ABSF (DC*EPS»BE TAD )-UN2)  130,140,140 

130  CALL  RANUKB(PP) 

JAB- J _ 

IF  (PP-P)  160,140,140 
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140  CALL  GRNUMBIALPHA)  _  _ _  _ _ 

CAU  6ft NUMB  {GAMMA  » 

_ RC -A»RC*S  I  GR*  SQft  TF  ( 1 . O-A—2  )  •ALPHA _ _ _ 

ISO  DC-6*DC*  S I GO* SQft TF 1 1 . 0-B*  *2 ) *GAMMA 

_ GO_TO_170 _ _ 

160  MM-NH41 

K-FLOATF  <  JAB » /DN*0 . 9999  _  _ _ 

JJ(K)-JJ(K)+1 

170  CONTINUE  _  _  _  _  _ 

FH-FLOATF INM) /F 
00  180  1-2, N 
Jjm-JJ(IUJJU-l) 

180  FFU.1-ELQATF  1JJUU/F _  _ _ 

FFIll-FLOATFI JJ(1» l/F 

_ PRINT  190*  A.S1GR. BETAR.B.S1GD.BETAD  _ 

190  FORMAT  (IHl,  49X,  11HC0RRELAT1 0N.6X,  9 HAIM  ERROR,  4X,  9HBALLISTIC/ 
133Xf  1 OH CQflftD I  MATE .  7X.  9HCDNSTANTS.  AX.  6HIM1LSI.  4X.  17HDISPERS1 
2 ON  (HILSI//34X,  5H RANGE*  13X,  F5.3,  10X*  F5.1,  9X.  F5.I/34X,  10H0E 

_ 3FLECTIQN*  BX*  1QX*  FS.l.  91U  Ei.ll _  _ 

PRINT  200,  C 

200  FORMAT  I/36X.  2SH  AIRCRAIT  SPEED  IKTAS1  -.F6.QI  _ 

PRINT  209,  R 

205  FORMAT  I34X.  2SH  FIRING  RATE  (RDS/M1N1 _ ■xfA.Q-L . . 

FAINT  210,  S 

210  FORMAT  134X.  43H  SLANT.  JL4NGE  AT  COMMENCEMENT  OF  FIRING  RUN  _ 

16HIFT )  >,F7.0I 

_ KM-FM . .  _ _  _ _ _ 

PRINT  215  ,  KN 

213  FORMAT  (34X.  25H  NO.  OF  ROUNDS  PER  PASS  -,161 _ _ _ 

PRINT  220,  P 

21Q  FJMMLJL14X,  3IH  CONDITIONAL  KILL  PROBABILITY  -,F6, 31 
PRINT  223,  FL 

225  FORMAT  </34X,  23H  TARGET  LENGTH  {FEET)  -,F7.l) 

PRINT  230,  M 

23Q  FORMAT  I34X,  23H  TARGET  KIQTH  (FEET) _ -tF7.il _ _ _ 

KF-F 

PRINT  235,  KF 

235  FORMAT  I/34X,  35M  NUMBER  OF  MONTE  CARLO  ITERATIONS  -,I61 

asNS-oyj)  _  . . 

PRINT  240,  NEMP 

240  FORMAT  I34X.  *3H  NUMBER  OF  EMPTY  PASS E1_T HRQUGH  RAjffiflMJHU _ 

1 12HGENER ATOMS  -.17) 

_  PRINT  249,  FH  ... _ _ _ _ _  ..  . . 

245  FORMAT  I///34X,  29H  PROBABILITY  OF  TARGET  KILL  »,F4.31 

TAINT  230  ...  .... 

250  FORNAT  I///5BX,  5H  KILL/34X,  13HN0.  OF  ROUNDS, 9X,  I IMPROBABILITY  I 

. . BiLlijLiiLtN _ _ _ _ _ 

ION-ON 

(R-IMON 

260  PRINT  2T0,  IR,  FF ( I > 

270  FORMAT  I39X,  15,  15X,  F5.3I 
GO  TO  10 
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APPENDIX  C 
DATA  SUBROUTINE 


1 .  Introduction: 

Many  computer  programs  require  the  flexibility  of  varying  any  or  all  of  the 
parameters  in  a  computer  run.  Although  FORTRAN  is  fairly  flexible  in  its 
arithmetic  and  control  statements,  its  input-output  statements  are  quite  rigid. 

In  order  to  read  cards  for  instance,  considerable  effort  must  be  expended  by  the 
FORTRAN  programmer  in  writing  his  input  statements.  This  subroutine  eliminates 
some  of  that  tedium.  The  concept  of  a  "data  set"  is  used.  A  data  set  consists  cf  a 
sequence  of  punched  cards  terminated  by  one  blank  card.  A  parameter  deck  for  a 
computer  run  may  consist  of  several  data  sets.  Such  a  parameter  deck  is  terminated 
by  two  blank  cards. 

2 .  Parameter  Addresses: 

The  primary  advantage  of  this  subroutine  over  FORTRAN  input  statements  re¬ 
sults  from  the  use  of  "parameter  addresses."  An  address  is  a  relative  location  in 
the  computer  memory.  It  is  the  subscript  of  an  array  or  matirx.  For  example,  in 
an  array  called  X  ,  the  parameter  value  X^  would  be  located  at  address  53  .  By 

using  the  parameter  addresses,  a  user  of  the  program  need  submit  only  those  param¬ 
eter  values  in  a  data  set  that  are  different  from  those  in  the  previous  set. 

Three  types  of  addresses  are  permitted  by  this  subroutine, 

(1)  A  numeric  address  consisting  of  one  to  five  characters,  each  of  which 
is  a  digit  0  -  9.  Such  an  address  (n)  refers  to  the  nth  element  in  a 
specified  array. 

(2)  An  alpha  address  consistir0  of  one  to  six  characters,  the  first  of  which 
must  be  alphabetic  (A-Z).  The  remaining  may  be  alphabetic  or  numeric 
(A-Z  or  0-9).  Such  an  address  refers  to  the  n™  element  in  a  specified 
array  (1  s  n  <  26),  where  the  first  character  of  the  address  corresponds 
to  n  as  the  26  letters  of  the  alphabet  correspond  to  the  integers  1-26. 

(31  A  matrix  address  consisting  of  two  or  more  numeric  fields  separated  by' 
commas.  For  example,  the  address  53,  47  refers  to  the  element  in  the 
53rd  row  and  the  47^  col  nn  of  a  two-dimensional  matrix.  There  is  no 
limit  to  the  number  of  dimensions  in  a  matrix  address. 

3.  Input  Card  Format: 

A  standard  submittal  form  (see  attachment)  has  been  designed  for  the  analyst. 
This  form  provides  for  entering  parameter  values  with  their  associated  addresses. 
The  user  indicates  blank  cards  to  separate  data  sets.  The  keypunch  operator  has 
the  option  of  punching  one  address  and  value  per  card,  or,  if  the  addresses  are 
sequential,  of  punching  one  address  and  several  values  on  a  card. 
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Only  columns  1-72  of  a  card  are  used.  Each  column  must  contain  one  of  the 
following:  a  digit  (o-9),  a  "+"  or  "  -  "  sign  or  a  dash,  a  letter  (A~Z),  a  period, 
a  comma,  or  a  blank.  Each  punched  card  must  contain  one  parameter  address 
The  address  may  start  in  column  1,  or,  if  desired,  may  start  in  a  later  column, 
provided  all  columns  before  it  are  blank.  The  address  is  terminated  by  at  least 
one  blank  column  Only  one  address  is  permitted  on  the  card.  Succeeding  col¬ 
umns  contain  one  or  more  parameter  values,  each  separated  by  one  or  more 
blank  columns  A  value  may  be  signed  or  unsigned.  The  length  of  the  value  field 
is  variable  No  blanks  are  permitted  within  a  value  field.  A  value  may  be  punched 
with  or  without  an  exponent.  An  exponent  is  recognized  by  the  presence  of  a  plus 
or  minus  sign  (or  dash)  between  the  fractional  part  and  exponent  part  of  the  value. 
Decimal  points  (periods)  may  be  punched  in  either  the  fractional  or  exponent  parts 
of  a  value.  If  more  than  one  value  is  punched  on  a  card,  those  after  the  first  will 
be  entered  at  sequential  addresses  relative  to  the  address  of  the  first  value. 

4  Usage: 

A  data  set  is  read  by  the  use  of  the  statement: 

CALL  DATA  (X,  I) 

in  a  FORTRAN  program  for  the  IBM  7090.  The  argument  X  is  the  name  of  an 
array  in  the  program.  The  argument  I  is  an  indicator  set  by  the  subroutine 
This  indicator  may  be  tested  by  the  main  program  upon  return  from  the  sub¬ 
routine  It  will  have  a  value  of  0  or  1  or  2. 

0:  The  subroutine  has  read  a  data  set.  The  main  program  will  normally 
proceed  to  operate  on  this  data. 

1:  The  subroutine  has  read  die  second  blank  card  which  terminates  the 
parameter  deck  The  main  program  will  normally  terminate  at  this 
point. 

2:  The  subroutine  has  read  a  "bad  '  data  card.  The  main  program  may 
terminate  the  run,  or  ignore  the  card  and  return  to  the  subroutine  to 
read  the  rest  of  the  data  set. 

If  the  cards  to  be  read  contain  matrix  addresses,  auJitional  arguments  must  be 
included  in  the  FORTRAN  calling  statement: 

CALL  DATA  (X,  D,,  D0,  D, . D  ,  I) 

izo  n 

where  Dj  is  the  i^  dimension  of  the  matrix  X. 

5  Method: 

See  the  attached  flow  chart.  DATA  reads  parameter  values  and  loading  addres¬ 
ses  from  cards.  If  sense  switch  5  is  up,  it  will  read  the  values  and  addresses  from 
tape  (unit  A2).  It  converts  the  values  to  floating  point  numbers,  and  stores  them  as 
elements  of  an  array  specified  in  the  calling  statement.  The  elements  are  specified 
by  the  addresses.  If  a  card  (or  tape  record)  Is  read  which  contains  non-permitted 
characters  (see  input  card  format  above),  DATA  prints  the  statement  "bad  data  card, " 
followed  by  an  image  of  the  card  itself. 
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6.  Coding  Information: 

See  die  symbolic  listing  included  in  this  appendix.  DATA  is  written  in  die  7090 
FAP  language.  It  must  be  used  in  conjunction  with  the  BELL  system.  It  requires 
401  words  storage  space. 
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SYMBOLIC  LISTING 


DATA 


A 1 
A 


A2 


8 


FAP 

entry  data 

SXA 

XI. 1 

SXA 

X2.2 

SXA 

X4.4 

CAL 

1.4 

ADD 

CORE 

STO 

XLOC 

AXT 

1.1 

SXA 

*+1.1 

CAL 

*#.4 

ANA 

MASK 

TNZ 

*+2 

TXI 

*-4.1.1 

SXA 

EX1T.1 

Txr 

*+1.1. ~1 

SXA 

*♦1.1 

CLA 

**.4 

STA 

A1 

STA 

F1A 

STA 

112 

AXT 

1.1 

SXD 

**.l 

TSX 

HHREA0.4 

PZE 

CARD 

TP* 

EXIT 

TRA 

BAD 

ST2 

ADORES 

ST2 

VALUE 

ST2 

ExP 

STZ 

POINT 

AXT 

1.1 

SXA 

FIELD. 1 

AXT 

13.1 

TNX 

HH.1.1 

AXT 

42.2 

SXA 

COLUMN. 2 

LXA 

COLUMN. 2 

TNX 

*2*2.6 

SXA 

COLUMN. 2 

LDO 

CARD* 12.1 

ROL 

36.2 

PXD 

0.0 

LGL 

6 

STO 

CHARAC 

ORA 

FLOAT 

FAD 

FLOAT 

STO 

NUMB 

AXT 

42.4 

CLA 

CHARAC 

CAS 

TABLE+42 .4 

TRA 

•+2 

TRA 

*♦3 

T I X 

*  3.4.1 

RETURN  INDICATOR  -  1 
READ  A  CARD 


ADDRESS  ■  0 
VALUE  -  0 
EXP  *  0 
POINT  ■  0 
FIELD  ■  1 


COLUMN  GT  72 


COLUMN  ■  COLUMN*! 
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TRA  PAD 

LXA  FIELD. 2 

TRA  Fl+1,2 

TRA  F7 

FIELD-7  (EXPONENT  FIELD) 

TRA  F6 

FIELD-6  (VALUE  FIELD) 

TRA  F5 

FIELD-5  (BLANKS  AFTER  ADDRESS) 

TRA  FA 

FIELD-4  (ARRAY  ADDRESS) 

TRA  F 3 

F J ELD-3  (ALPHA  ADDRESS) 

TRA  F2 

FIELD-2  (NUMERIC  ADDRESS) 

FI 

TXH  B.A.A1 

FIELD-1  (BLANKS  BEFORE  ADDRESS) 

F1A 

STZ  ** 

RETURN  INDICATOR  •  0 

AXT  2,2 

FIELD  -  2 

SKA  FIELD. 2 

TXH  H.4.31 

NUMERIC  CHARACTER 

TXH  1,4.28 

SIGN  OR  DASH 

TXL  ttAD.4.2 

AXT  3.2 

ALPHA  CHARACTER.  FIELD  -  3 

SXA  FIELD. 2 

TXI  **1 .4,-2 

ADDRESS  -  NTH  ALPHA 

F2 

SXA  ADORES. A 
TRA  B 

TXH  J.A.Al 

BLANK  CHARACTER 

TXH  H.A.31 

NUMERIC  CHARACTER 

TXH  K.4.28 

SIGN  OR  DASH 

TXH  BAD. A. 2 

TXH  L.4.1 

PERIOD 

AXT  4.2 

COMMA. FIELD  •  4 

SXA  FIELD. 2 

AXT  2,2 

SXA  DIMENS, 2 

DIMENSION  -  2 

cla  ADORES 

TEST  ADDRESS 

F2A 

TZE  BAD 

TM1  BAD 

STZ  ADDER 

ADDER -C 

F  3 

TRA  P 

TXH  J.4,41 

BLANK  CHARACTER 

TXH  B.4,31 

NUMERIC  CHARACTER 

TXH  K.4.28 

SIGN  OR  DASH 

TXH  B.4.2 

ALPHA  CHARACTER 

TXH  L.A.l 

PERIOD 

FA 

TRA  BAD 

TXH  M.4,41 

BLANK  CHARACTER 

TXH  N.4,31 

NUMERIC  CHARACTER 

TXH  P.A.28 

SIGN  OR  DASH 

TXH  BAD  » A  .2 

TXH  <3,4.1 

PERIOD 

TRA  T 

COMMA 

F5 

TXH  B.4,41 

BLANK  CHARACTER 

AXT  6.2 

FIELD  •  6 

SXA  FIELD. 2 

TXH  U.4.31 

NUMERIC  CHARACTER 

TXH  M.4,28 

SIGN  OR  DASH 

TXH  BAD, 4, 2 

TXH  G.4,1 

PERIOD 
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TRA  BAO 

F6 

TXH  X,4. .1 

BLANK  CHARACTER 

TXH  U.4.31 

NUMERIC  CHARACTER 

TXH  7,4,2# 

TXH  BAD, 4. 2 

SIGN  OR  DASH 

TXH  G.t.l 

TRA  BAD 

PERIOD 

F7 

TXH  2,4,41 

BLANK  CHARACTER 

TXH  BB.4,31 

NUMERIC  CHARACTER 

TXH  EE, 4, 28 

TXH  BAD, 4. 2 

TXt.  SAD. 4,1 

SIGN  OR  DASH 

6 

AXT  1.2 

SXA  POIHT.2 

TRA  8 

PERIOD.  POINT  •  1 

H 

LOO  ADORES 

MPT  HlO 

XCA 

ACL  CHARAC 

STO  ADORES 

TRA  B 

ADDRESS  ■  10  X  ADDRESS 

I 

TXH  6.4.30 

♦  SIGN 

CLA  ADORES 

SSH 

STO  ADORES 

TRA  B 

SET  SIGN  Or  ADDRESS 

J 

AXT  5.2 

SXA  FIELD. 2 

TRA  B 

PIELD  •  5 

K 

TXH  LI. 4. 30 

*  SIGN 

CLA  VALUE 

SSH 

STO  VALUE 

TRA  LI 

SET  SIGN  OF  VALUE 

L 

AXT  1,2 

SXA  POINT. 2 

POINT  •  1 

LI 

AXT  A .2 

SXA  FIELD. 2 

TRA  e 

FIELD  •  6 

H 

AXT  5.2 

SXA  FIEL0.2 

TRA  S 

FIELD  »  5 

N 

LOO  ADDER 

NPY  HlO 

STO  ADDER 

TSX  T1.4 

HPT  CHARAC 

XCA 

ADO  ADDER 

STO  AOOER 

TRA  S 

ADOER  -  10  X  ADDER  ♦  N 

P 

TXH  R.4.30 

♦  SIGN 

C<  A  VALUE 

SSH 

SET  SIGN  OF  VALUE 
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STO  VALUE 
TRA  R 

Q  AXT  1,2  POINT  ■  1 

SXA  POINT *2 

R  AXT  6,2  FIELD  ■  6 

SXA  FIELD, 2 

S  LXA  EXIT, 2  CHECK  DIMENSION 

TXI  1 *2  »-3 
PXA  0,2 
SUB  DIMENS 
TNZ  8AD 

T  TSX  Tl,4  ADOER* ADDER -PROD 

Cl  A  ADDER 

SUB  PROD 
STO  ADDER 

TZE  BAD  CHECK  ADDER 

TMI  BAD 

ADD  AODRES 

STO  AODRES 

CLA  DIMENS 

ADD  HI 

STO  OIMENS 

TRA  F2A 

T1  SXA  T4»4  PROD  ■  PRODUCT  OF  DIMENSIONS 

CLA  HI 
STO  PROD 
STA  T3 

LXA  OIMENS, 2 
TXI  *41 ,2,-1 
LXA  X4.4 
T 2  CAL  TS 

ADO  HI 
STA  T3 

T 3  CLA  *«,4 

STA  *-,1 
LDQ  ** 

ROL  It 
MPV  PROD 
STO  PROD 
T  I X  T  2 , 2  «  1 

T*  AXT 

TRA  1,4 

U  CLA  POINT  TEST  POINT 

TNZ  V 

LD3  VALUE  VALUE  •  10  X  VALUE  ♦  N 

FMP  DEC 10 

SSP 

FAD  NUMB 
LOO  VALUE 
LLS  0 
STO  VALUE 

TRA  B 

V  LXA  POINT ,4  VALUE  •  VALUE  ♦  N/ i iO**POINT ) 

CLA  NUMB 
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FDP 

DEC  10 

XCA 

TIX 

•-2.4.1 

LDQ 

VALUE 

LLS 

0 

FAD 

VALUE 

STO 

VALUE 

TRA 

DO 

w 

TXH 

B.4.30 

CLA 

VALUE 

SSM 

STO 

VALUE 

TRA 

B 

X 

CLA 

XLOC 

SUB 

ADDRES 

STA 

*♦2 

CLA 

VALUE 

STO 

M 

TRA 

AA 

Y 

AXT 

7.2 

SXA 

HtLD.2 

TXH 

6G.4.30 

CLA 

EXP 

SSM 

STO 

EXP 

TRA 

GG 

z 

CLA 

XLOC 

SUB 

ADDRES 

STA 

Z1 

CLA 

DEC  10 

LOO 

EXP 

CALL  EXPO 

XCA 

FWP 

VALUE 

Z1 

STO 

## 

AA 

STZ 

VALUE 

TRA 

FF 

08 

CLA 

POINT 

TNZ 

CC 

LDQ 

EXP 

FMP 

DEC10 

S$P 

FAO 

NUMB 

LOQ 

EX® 

LLS 

0 

STO 

EXP 

TRA 

B 

CC 

LXA 

POINT  *4 

CLA 

NUMB 

FOR 

DEC10 

XCA 

TIX 

*-2.4*1 

LDQ 

EXP 

LLS 

0 

*  SIGN 

SET  SIGN  OF  VALUE 


X( ADDRESS)  ■  VALUE 


FIELD  •  2 
+  SIGN 

SET  SIGN  OF  EXP 


X( ADDRESS)  ■  VALUE  X  10**£XP 


value  ■  0 

TEST  POINT 

EXP  •  10  X  EXP  ♦  N 


EXP  •  EXP  *  N/)10**POINT> 
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DO 

FAu  EXP 

STO  EXP 

CL*  POINT 

POINT  ■  POINT  ♦  1 

EE 

ADO  HI 

STO  POINT 

TRA  8 

CLA  XLOC 

X ( ADDRESS )  -  VALUE  X  10«»EXP 

EE1 

SUS  ADORES 

STA  EE1 

CIA  DEC 10 

LDO  EXP 

CALL  EXPO 

XCA 

FMP  VALUE 

STO  *• 

PXD  0.0 

VALUE  ■  0 

TXH  •♦2.4,30 

♦  SIGN 

SSN 

SET  SIGN  OF  VALUE 

FF 

STO  VALUE 

STZ  EXP 

EXP  •  0 

AXT  5.2 

FIELD  ■  5 

SXA  FIELD, 2 

CAL  AODRES 

ADDRESS  ■  ADDRESS  ♦  1 

GG 

ADO  HI 

SLW  AODRES 

STZ  POINT 

POINT  •  0 

HH 

TRA  6 

LXA  FIELD. 1 

TXL  JJ.1.1 

FIELD-1,  EXIT 

TXL  BAD, 1,4 

TXL  A, 1,5 

FIELD-5.  READ  ANOTHER  CARD 

TXH  11.1,6 

CLA  XLOC 

FIELD-6.  XI  ADDRESS)  •  VALUE 

II 

SUB  AODRES 

STA  **2 

CLA  VALUE 

STO  •* 

TRA  A 

TXH  BAD, 1,7 

CLA  XLOC 

FIELD-7, 

SUB  AODRES 

XI  ADDRESS)  •  VALUE  X  10*»EXP 

III 

STA  III 

CLA  0EC10 

LOO  EXP 

CALL  EXPO 

XCA 

F*P  VALUE 

STO  ** 

8*0 

TRA  A 

TSX  HPRINT.4 

112 

PZE  PRINT, 0,15 
AXT  2,1 

SXD  •*» 1 

XI 

AXT  ••»! 

C- 12 


INTERIM  RESEARCH  MEMORANDUM 

IRM-43 


X2  AXT  *»*2 
X4  AXT  «MI»4 
EXIT  TRA  **.* 

MASK  OCT  777777700000 
PRINT  SCO  3  BAD  DATA  CARO... 
CARD  BSS  12 
ADORES  HiR  ** 

VALUE  HTR  •• 

EXP  HTR  *• 

POINT  HTR  ** 

FIELD  HTR  ** 

COLUMN  HTR  ** 


TABLE  OCT  60  BLANK 

OCT  0  0 

OCT  1  1 

OCT  2  2 

OCT  3  3 

OCT  *  * 

OCT  5  3 

OCT  6  6 

OCT  7  7 

OCT  10  8 

OCT  11  9 

OCT  20  ►  SI6N 

OCT  40  -  SIGN 

OCT  14  DASH 

OCT  71  Z 

OCT  70  Y 

OCT  67  X 

OCT  66  ¥ 

OCT  65  V 

OCT  64  U 

OCT  63  T 

OCT  62  5 

OCT  51  R 

OCT  30  0 

OCT  47  P 

OCT  46  0 

OCT  45  N 

OCT  44  M 

OCT  43  L 

OCT  42  K 

OCT  41  J 

OCT  31  1 

OCT  30  H 

OCT  27  G 

OCT  26  F 

OCT  25  E 

OCT  24  D 

OCT  23  C 

OCT  22  B 

OCT  21  A 

OCT  S3  PERIOD 

OCT  73  COMMA 
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CHARAC  HTR  *♦ 

DIMENS  HTR  ** 
aOOER  HTR  ** 

HIO  HTR  iO 

DECIO  DEC  10.0 

HI  HTR  1 

PROD  HTR 

AMASK  OCT  77777 

FLOAT  OCT  233000000000 

NUMB  HTR  *• 

XLOC  HTR  ** 

CORE  OCT  100001 
JJ  SYN  XI 
END 
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OEG  COMPUTER  DATA  SUBMITTAL  FORM 


Address  Value 


NOTES: 


A  value  of  zero  must  be  entered  as  0,  not  left  blank. 

Decimal  pts.  may  be  omitted  If  understood  to  follow  the  rightmost  digit. 

The  value  3  X  10  5  may  be  entered  as  .00003  or  3-S,  not  as  3  X  10~‘\ 
The  factor  portion  of  a  value  may  not  contain  more  than  8  digits. 

The  exponent  portion  of  a  value  must  lie  within  the  range  ±39. 

Exponents  may  be  omitted  if  zero.  If  not,  they  must  be  signed. 

Blank  cards  should  be  indicated  by.  - b - . 
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APPENDIX  D 
GRNUMB  SUBROUTINE 


1.  Purpose: 

GRNUMB  provides  a  floating  point  pseudo-random  number  X.  The  distribution 
of  successive  values  of  X  are  Gaussian  with  a  mean  of  zero  and  a  standard  devia¬ 
tion  of  one. 


2.  Method: 

Consider  the  set  of  uniformly  distributed  pseudo -random  numbers  .  GRNUMB 

generates  a  sequence  of  Y.  by  the  method  of  congruences: 


Y.  =  2'35’(515235Y._1  ,  mod  235) 
over  the  range  0  i  Y .  <  1  .  The  variance  of  this  uniform  set  is 

ay2  =  J1(Y-l/2)2dY  =  1/12  . 

o 

If  X  is  the  mean  of  any  selection  of  m  of  the  uniform  numbers  Y  ,  the  Central 
Limit  Theorem  states  that  the  variable  X  approaches  a  normal  distribution  where 
m  is  sufficiently  large.  A  satisfactory  value  for  m  is  30.  Values  of  X  are 
generated  as  a  sequence  of  XR  ,  where  n  denotes  the  n^1  entry  to  GRNUMB. 

- -  m  30 

l  <vf  1/2>  -nl  <V1/2> 

i=l  i=l 

where  Y  =  X  .  ,  and  X  =  2-^  .  The  variance  of  this  normal  set  is  1 . 
on-.'  o 

3.  Usage: 

X  is  obtained  by  use  of  the  statement: 


•K 

% 

i 

I 

* 

I 

I 


CALL  GRNUMB  (X) 


in  a  FORTRAN  program  for  the  IBM  7090. 

4.  Coding  Information: 

See  the  symbolic  listing  on  the  following  page.  GRNUMB  is  written  in  the 
7090  FAP  language.  It  requires  40  words  storage  space  and  900  microseconds 
operating  time. 
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SYMBOLIC  LISTING 


GRNUMB 


C 


D 

E 


F 

XXI 

NUM 

NUMB 

MULT 

CHAR 

MAGIC 

H6 

H125 


FAP 

REM  GRNUMB  G.  WESTLUND  18  JUNE  1962  (7090) 

REM  GAUSSIAN  DISTRIBUTED  RANDOM  NUMBER  GENERATOR. 

REM  ENTER  VIA  FORTRAN  STATEMENT  CALL  GRNUMB  (X) 

REM  SEQUENCE  STARTS  AT  DEC  1,  YIELDS  X  WITH  STAND.  DEV.  =1 . 

ENTRY  GRNUMB 

SXA  XXI,  1 

CLA  1,  4 

STA  F 

AXT30,  1 

STZNUM 

LDQ  NUMB 

MPY  MULT 

STQ  NUMB 

CLA  NUMB 

SUB  CHAR 

ARS  4 

ADD  NUM 

STO  NUM 

TIXC,  1,  1 

LDQ  NUM 

MFY  MAGIC 

LRS  27 

TZE  D 

LRS  8 

CLA  HI 25 

APDH8 

LLS  8 

ALS  19 

TRA  E 

CLA  H125 

ALS  27 

STO  NUM 

CLA  HUS 

LLS  27 

FAD  NUM 

STO  •• 

AXT  **,  1 
TRA  2,  4 
HTR  ♦* 

DEC  1 

DEC  30517578125 

TDC  0,  0.  0 

DEC  0. 31622780B0 

DEC  8 

DEC  125 

END 
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APPENDIX  E 
RANUMB  SUBROUTINE 


1.  Purpose: 

RANUMB  provides  a  floating  point  pseudo -random  number  X.  Successive 
values  of  X  are  uniformly  distributed  over  the  range  0  s  X  <  1 . 

2.  Method: 

The  generated  value  of  X  is  a  member  of  a  set  of  pseudo-random  numbers. 
This  set  is  a  sequence  erf  Xr  ,  where  n  denotes  the  n^1  entry  to  RANUMB.  The 

set  is  generated  by  the  method  of  congruences: 

Nn  =  (S15Nn_i),  mod235 

\  35 

\  X  =  N  /2" 

'■  n  n 

where  N  =  1 . 
o 

3.  Usage: 

X  is  obtained  by  use  of  the  statement: 

CALL  RANUMB  (X) 

in  a  FORTRAN  program  for  the  IBM  7090. 

4.  Coding  Information: 

See  the  symbolic  listing  on  the  following  page.  RANUMB  is  written  in  the 
7090  FAP  language.  It  requires  15  words  storage  space  and  41  microseconds 
operating  time. 
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SYMBOLIC  LISTING 

FAP 

REM  RANUMB  G.  WESTLUND  18  JUNE  1962  <7(990) 

REM  UNIFORMLY  DISTRIBUTED  RANDOM  NUMBER  GENERATOR. 
REM  ENTER  VIA  FORTRAN  STATEMENT  CALL  RANUMB  <X) 
REM  SEQUENCE  STARTS  AT  DEC  1.  YIELDS  (0  LE  X  LT  1). 
ENTRY  RANUMB 
RANUMB  CLA  1.  4 
CTA  r 

LDQ  NUMB 
MPY  MULT 
STQ  NUMB 
PXD0,  9 
RQL  9 
LGL  27 
ACL  CHAR 
FAD  CHAR 
C  STO  ** 

TRA  2.  4 

NUMB  DEC  1 

MULT  DEC  3*1517578125 

CHAR  TIX  0,  0,  0 

END 


E-2 
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APPENDIX  F 

SAMPLE  PROBLEM  SUBMITTAL  FORM 


OEG  COMPUTER  DATA  SUBMITTAL  FORM 

Submitted  by.  John  j)p£ _  Date:  Z(p  Jllltj  I  9 (p 3 

Program  No. _ Eat.  Time  ^  H»in _  Classification  ^yjcJfljjSjjrgj 


Special  Instructions: 


Address 

Value 

Address 

Value 

Address 

Value 

Address 

Value 

1 

0.8 

£ 

0.8 

3 

b 

4 

4 

5 

2 

b 

1 

7 

5000 

& 

1000 

9 

500 

10 

100 

1 1 

0.2 

12 

10 

13 

10 

14 

2000 

15 

0 

1  Co 

10 

-b 

-b . 

NOTES: 

1.  A  value  of  zero  must  be  entered  as  0,  not  left  blank. 

2.  Decimal  pts.  may  be  omitted  if  understood  to  follow  the  rightmost  digit. 

3.  The  value  3  X  10  5  may  be  entered  as  .  00003  or  3-5,  not  as  3  X  10"5. 

4.  The  factor  portion  of  a  value  may  not  contain  more  than  8  digits. 

5.  The  exponent  portion  of  a  value  must  lie  within  the  range  ±39. 

6.  Exponents  may  be  omitted  if  zero.  If  not,  they  must  be  signed. 

7.  Blank  cards  should  be  indicated  by  - b - . 


F  - 1 

(REVERSE  BLANK) 
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APPENDIX  G 

SAMPLE  PROBLEM  OUTPUT 


CORRELATION 

COORDINATE  CONSTANTS 


AIM  ERROR  BALLISTIC 
(MILS!  DISPERSION  (MILS) 


RANGE  0.800 

OEFLECTI ON  0.800 


£.0  2.0 

4.0  1.0 


AIRCRAFT  SPEED  (KTAS)  *  500. 

FIRING  RATE  (RDS/MIN*  »  1000. 

SLANT  RANGE  AT  COMMENCEMENT  OF  FIRING  RUN  (FT)  «  5000. 

NO.  OF  ROUNDS  PER  PASS  *  100 

CONDITIONAL  KILL  PROBABILITY  -  0.200 

TARGET  LENGTH  (FEET)  *  10.0 

TARGET  WIDTH  (FEET)  -  10.0 

NUMBER  OF  MONTE  CARLO  ITERATIONS  *  2000 

NUMBER  OF  EMPTY  PASSES  THROUGH  RANDOM  NO.  GENERATORS  -  0 


PROBABILITY  OF  TARGET  KILL  *  0.993 


NO.  OF  ROUNDS 
10 
20 
30 
40 
50 
60 
70 
80 
90 
100 


KILL 

PROBABILITY 

0.048 

0.114 

0.200 

0.287 

0.384 

0,499 

0.658 

0.812 

0.950 

0.993 


G-l 

(REVERSE  BLANK) 


